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Abstract 

We investigate the potential to observe double parton scattering at the Large Hadron Collider 



^1^' in pp — )• WbbX — )• iubbX at 7 TeV. Our analysis tests the efficacy of several kinematic variables 

j-^. in isolating the double parton process of interest from the single parton process and relevant 

^vq ! backgrounds for the first 10 fb^^ of integrated luminosity. These variables are constructed to 

(^ ' expose the independent nature of the two subprocesses in double parton scattering, pp — t- ivX and 

vn ■ 

pp —7- bbX. We use next-to-leading order perturbative predictions for the double parton and single 

cn 

r^ . parton scattering components of Wbb and for the pertinent backgrounds. The next-to-leading order 

o ■ 

contributions are important for a proper description of some of the observables we compute. We 
find that the double parton process can be identified and measured with significance S/yB ~ 10, 
S^ . provided the double parton scattering effective cross section a^s ~ 12 mb. 
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I. INTRODUCTION 

The successful operation of the Large Hadron Colhder (LHC) and its detectors opens a 
new era in particle physics. The higher energies and larger luminosities at the LHC make it 
possible to explore new physics scenarios and to investigate unexplored aspects of established 
theories such as quantum chromo dynamics (QCD). 

The standard picture of hadron-hadron collisions is shown on the left side of Fig. [1] One 
parton from each proton partakes in the hard scattering to produce the final state. The 
probability density for finding parton i in a proton with momentum fraction Xi and at the 
factorization scale /i is parametrized by the parton distribution function (PDF) fHxi,^). 
In this single parton scattering (SPS) scenario, the differential hadronic cross section neatly 
factors into: 

d(^PP^ " 5Z / /p(^i' /^)/pK' lAdc^iA^ii x[,fi)dxidx[ . (1) 

The "short-distance" partonic cross section daij is computed in perturbation theory, whereas 
the PDFs are nonperturbative objects and must be extracted from experiment. 

This simple picture of proton-proton collisions is incomplete. The full description of 
hadronic collisions involves other elements including initial- and final-state soft radiation, 
underlying events, and multiparton interactions. Double parton scattering (DPS) describes 
the case in which two short-distance subprocesses occur in a given hadronic interaction, 
with two initial partons being active from each of the incident protons. The general picture 
of DPS is shown on the right side of Fig. [H Given the small probability for single parton 
scattering in hadronic collisions, it is often assumed that the effects of double (or multiple) 
parton scattering may be ignored or subsumed into the parametrization of the underlying 
event. Nevertheless, it is worth exploring theoretically and investigating experimentally 
whether a second distinct hard component may be identified in events at the LHC. Some 
evidence for DPS has been observed at the CERN Intersecting Storage Rings Jl| 
Super Proton Synchrotron [2i], and more recently, at the Fermilab Tevatron 

In an earlier study Sj, we investigated the DPS and SPS contributions at the LHC 
to the four-parton final state pp — )■ bbjjX in which a bb system is produced along with 
two jets j. We showed that there are characteristic regions of phase space in which the 
DPS events are expected to concentrate, and we developed a methodology to measure the 
effective size of DPS. Precise measurements of DPS at the LHC will provide insight into 






FIG. 1: Schematic depiction of single parton scattering (left) and double parton scattering (right). 
In single parton scattering, one parton from each hadron is active in the scattering and the partonic 
process is A{ij — )• abed). Double parton scattering assumes two partons from each hadron are active 
in the hard scattering and the total partonic process consists of two independent subprocesses 
A{ij — ;> ab) and A{k£ — )- cd). 

parton correlations, nonperturbative dynamics in hadron-hadron collisions, the structure of 
the proton, and parton distribution functions. 



Theoretically, the study of DPS phenomenology has a long history |5l-l37j]. Under the 
assumption of weak dynamic and kinematic correlations between the two hard-scattering 
subprocesses, the general approach in these studies is to assume the differential hadronic 
cross section takes a factored form in analogy to Eq. ([T]): 
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where m is a symmetry factor which is equal to 1 (2) if the two hard-scattering subprocesses 
are identical (nonidentical) . The joint probabilities Hp'''{xi,X2, fiA, /^b) can be approximated 
as the product of two single PDFs: 



H'/{Xi,X2,fiA,IJ'B) = rJxi,fiA)f^{x2,IJ,B) 



(3) 



Given that one hard scattering has taken place, the parameter aes measures the size of 
the partonic core in which the flux of accompanying short-distance partons is confined. 
Typical values in phenomenological studies focus on the 10-12 mb region, consistent with 
measurements from the Tevatron collider |3|, |j]. In writing Eqs. ([2]) and ([2]), we ignore 
possible strong correlations in longitudinal momentum. However, for the small values of x 
expected at the LHC, this should be a good approximation 



In order to observe a DPS signal, it is advantageous if: 1) the cross sections for the two 
individual processes which make up the DPS process are large (in the nib - pb range), and 2) 
the final state contains objects which can be easily tagged (or identified). The bbjj channel 
possesses both of these qualities. By focusing on distributions which contain kinematic 
information about the entire final state, we were able to show [5] that the DPS component 
of the final state can be observed despite the presence of a dominant SPS component over 
a wide kinematic range. 

In this paper, we examine another final state which may be a good candidate to observe 
DPS, namely the production oi a.W boson in association with a pair of bottom quark jets. 
In the DPS contribution to this final state, one hard scattering produces the W via the 
Drell-Yan mechanism, while the other hard scattering produces a bb pair. The cross sections 
for these two processes are individually large, and the charged lepton from the W decay 
(along with the bottom quarks in the final state) provides a relatively clean signal to tag 
on. Our purpose is to establish whether double parton scattering can be observed as a 
discernible physics process in Wbb production at LHC energies. We remark, however, that 
the Wbb final state is a significant background for the production of a Higgs boson H in the 



HW mode, where the Higgs boson decays into a pair of bottom quarks 2l| , and that it 



can be a background in channels where new physics may arise such as in single top quark 



production [38[. Once DPS production of Wbb is observed, it would be interesting to assess 
its potential significance as a background in such searches. A realistic study would require 
knowledge of the effective cross section for double parton scattering in the Wbb channel and 
a set of optimized physics cuts pertinent for the search in question. We leave this topic for 
possible future work. 

The rest of the paper is structured as follows. In Sec. [TTl we outline our procedure 
for computing the DPS and SPS contributions to Wbb at the LHC, and we discuss and 
evaluate backgrounds to the same final state. Our SPS and DPS event rates are computed 
at next-to-leading order with the aid of the POWHEG BOX code i39|. This section includes 
a specification of the basic acceptance cuts we use in defining the event sample. Section IIIII 
is devoted to the role of the large ti background. We find that a cut to eliminate events 
with large missing transverse energy is effective in suppressing this background. In Sec. 
llVt we focus on discrimination between the DPS and SPS contributions to Wbb. We study 
various single variable and two-dimensional kinematic distributions to bring out the DPS 



contribution more cleanly. By utilizing cuts that enhance the DPS Wbb sample, we find that 
the DPS signal can be observed with a statistical significance in the range S/yB ~ 12 — 15. 
Section |V] contains our summary. 

II. DPS AND SPS CONTRIBUTIONS TO Whh PRODUCTION AT THE LHC 

We begin with the premise that there are DPS and SPS components of the same Whb 
final state. Our aim is to try to pick out the DPS component and to study its distinct 
properties. In this section, we outline our method for computing event rates for Whb from 
DPS and SPS as well as the backgrounds for the same final state at the LHC. We perform all 
calculations at a center-of-mass energy of ^/s = 7 TeV. Event rates are quoted for 10 fb^^ of 
integrated luminosity. For the DPS case, Wbb production is computed using Eq. ([2]) where 
it is assumed that one hard scattering produces the W boson via the Drell-Yan mechanism 
{qq — > W"^ at leading order), while the other scattering produces the bb system (with either 
gg — )■ bb or qq -^ bb). Schematically, we can represent the partonic DPS process as: 

{ij -> W^) ® [kl -^ bb) . (4) 

The individual SPS processes which make up the DPS process in Eq. (j4]) are generated using 



the POWHEG BOX event generator 39l-l41l| which includes next-to-leading order (NLO) 



QCD corrections for both, plus shower emission. The ® symbol denotes the combination 
of one event from each of the W^ and hb final states. All events are produced using the 



two- loop evaluation of ^^(yu) (where /i = M^y for the W process and /i = \Jml + p^ for the 
bb production) and CTIO NLO PDFs |42|. 



In the SPS production of Wbb, one hard scattering produces the complete final state. The 



so generated using the POWHEG BOX [43] which implements 



events from this process are a.. 
the NLO calculation of Ref. J44|. 

Extracting evidence for DPS Whb production is complicated by the fact that many stan- 
dard model processes imitate the Wbb — )• bbiv final state. In particular, we consider contri- 
butions from the following final sates: 

• Top quark pair production tt where either {i) both t's decay semileptonically (denoted 
by t^t^), and one of the charged leptons is missed, or {ii) where one t decays semilep- 
tonically while the other decays hadronically (denoted by tith) and two jets are either 



missing or do not pass the threshold and isolation cuts. 

• Single top quark production {tb, ib, tj and ij modes) where t — )■ W~^b (t — )■ W~b). 

• Wjj, where both light jets are mistagged as a 6 jets. 

• Wbj where the light jet is mistagged as a 6 jet. 

We also considered the following processes, which have a negligible contribution after cuts: 

• bbj where one b quark gives an isolated lepton and the light jet is tagged as a 6 jet. 

• Zbb where one lepton from the Z decay goes missing. 

• bbbb {bbcc) production where at least one heavy quark gives an isolated lepton and the 
other does not pass the threshold cuts. 

The Wjj background (where both jets fake bottom quark jets) can be produced in both SPS 
and DPS processes. We compute the DPS contribution using the same method discussed 

n 

above for Wbb, using the additional POWHEG code [45], and we include it in our analysis. 



The top pair J41[ and single top [46|| SPS processes are also gene rated using the POWHEG 



BOX. Othe. SPS p.oce.es a.e generated u.ng MadEvent W o. ALPGEN M; W,J i. 

reweighted using a K factor obtained with MCFM |49| . 

In order to avoid soft and collinear divergences in the processes that we compute at LO, 
we apply a minimal set of generator-level cuts: 

PTj > 15GeV , \r],\ < 4.8 , \r]e\ < 2.5, (5) 

PT,b> 15GeV , \rib\ < 2.5 (6) 

ARj^b)j{b) > 0.4 , ARj^b)e > 0.4 , (7) 

where r] is the pseudorapidity and ARik is the separation in the azimuthal-pseudorapidity 
plane between the two objects / and k: 



AR, 



Ik 



\/{vi - mf + i<Pi - <Pkf . (8) 



Some of these generator-level cuts cannot be applied for processes computed with POWHEG, 
but they are applied subsequently to ensure equal treatment of all event samples. 



A. Simulation 

We concentrate on the final state in which there are two b jets, a hard lepton, and missing 
transverse energy ^y. To identify the Wbb final state and reduce backgrounds, we begin 
with simple identification cuts on the generated event samples. First, we consider only 
leptonic decays of the W boson (W — )■ iu). We focus on the case i = fi, since electrons with 
low transverse momentum can be easily faked by light jets. We limit the hadronic activity 
in our events to include exactly two hard jets, both of which must be identified as bottom 
quark jets. Finally, all events (DPS and SPS Wbb as well as backgrounds) are required to 
pass the following acceptance cuts: 

PT,6>20GeV, |r/fe|<2.5, (9) 

20GeV<pT,M <50GeV, Ir^^l <2.1, (10) 

^T > 20 GeV (11) 

AR,i > 0.4 , ARb^ > 0.4 . (12) 

The cut on the missing transverse energy ^j- > 20 GeV is motivated by the fact that 
the neutrino momentum in W decay is not observed. The 20 GeV cut on the b jets and the 
lepton is invoked to eliminate contributions from the underlying event. The upper lepton 
Pt cut is used to reject boosted W bosons, as in the case where a W boson originates from 
a t-quark decay, or when the W recoils against the bb pair in SPS. 

To account for b jet tagging efficiencies, we assume a 6-tagging rate of 60% for b quarks 
with pTfi > 20 GeV and \rih\ < 2.5. We apply the muon identification efficiencies found in 



the ATLAS Technical Design Report 50|- Detector resolution effects are accounted for by 
smearing the final-state energy according to: 

SE a , ,,„, 

where a = 50% and b = 3% for jets and a = 10% and b = 0.7% for leptons. Light jets (jets 
from u, d, s and c quarks as well as gluons) can "fake" bottom quark jets and we account 
for this by applying a mistagging rate for the gluon and u, d and s quarks of: 

euAs,g-.b = 0.67% (14) 



for pT,j < 100 GeV and: 

^u,d,s,g^b = 2% (15) 

for ptj > 250 GeV. For pt values between 100 and 250 GeV, we linearly interpolate the 
fake rates. Finally, for c quarks, we assume a fake rate of: 

ec^b = 10% (16) 

for pt^c > 50 GeV and we linearly interpolate fake rates for pj- < 50 GeV. 

Table [T] shows the number of events from the Wbb final state (DPS and SPS) and the 
backgrounds both before (column labeled "generator-level cuts") and after the acceptance 
cuts, detector effects, and mistagging effects are applied (column labeled "acceptance cuts"). 
In these results and those that follow, we sum the W~^ and W~ events. In evaluating the 
DPS processes, we assume a value o"eg ~ 12 mb for the effective cross section. However, we 
stress that the goal is to motivate an empirical determination of its value at LHC energies. 
The acceptance cuts are very effective against the Wjj final states, both for DPS and SPS. 
The results in Table |T] make it apparent that Wbb production from SPS and the top quark 
pair background are the most formidable obstacles in extracting a DPS signal. We address 
background rejection in the next two sections. 



III. tt BACKGROUND REJECTION 

We examine three possibilities to reduce the ti background: a cut to restrict J^t from 
above, rejection of events in which a top quark mass can be reconstructed, and a cut to 
restrict the transverse momentum of the leading jet. In the end, an upper cut on ^t in the 
event appears to offer the best advantage. Indeed, one would expect that ^t in Wbb events 
would be smaller than ^t in ti events. Top quark decays give rise to boosted VT^'s which, 
after decay, should result in larger values of missing Et compared to the Wbb process. The 
^T distribution is shown in Fig. [2] for the DPS component of Wbb, the SPS component of 
Wbb, and all backgrounds (left). On the right, we show the DPS component of Wbb and 
the ti background alone. The plot on the right shows that the DPS signal is produced in 
the region of relatively small ^t and the ti background has a harder spectrum in !Pt- One 
way to suppress the ti background while leaving the DPS signal unaffected is to impose 
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Process 


Generator-level cuts 


Acceptance cuts 


^T < 45 GeV 


S',^ < 0.2 


WHb (DPS) 


10000 


247 


231 


173 


W'^bb (SPS) 


44000 


1142 


569 


114 


tt 


225000 


1428 


290 


13 


W^jj (DPS) 


476000 


43.5 


37.7 


27.3 


W^jj (SPS) 


20300000 


101 


55.7 


19.6 


Single top 


20000 


492 


168 


15 


W'^bj 


153000 


152 


53.1 


8.2 



TABLE I: Numbers of events before and after the various cuts are applied for 10 fb^^ of data. After 
acceptance cuts, SPS Wbb production and tt production dominate the event rate. A maximum ^^ 
cut reduces the background from tt significantly. A maximum cut on S' improves the DPS/SPS 
ratio in Wbb production. 



a maximum ^t cut in the 40-60 GeV range. In the analysis that follows, we include a 
maximum ^^ cut of 45 GeV in addition to the acceptance cuts outlined above. 

The effects of the maximum ^t cut are shown in the fourth column of Table [11 This cut 
eliminates about 80% of the ti background that remained after the initial acceptance cuts. 
The cut is also effective at reducing the single top quark and Wbj backgrounds, eliminating 
about 67% in both cases. On the other hand, 93% of the DPS Wbb events and 46% of the 
SPS Wbb events are retained. 

Backgrounds from events that contain a real top quark, such as the ti and single top 
events, might be reduced if one could reconstruct a top quark mass distribution M^i, from 
the final-state objects (bottom quarks, charged leptons, and neutrinos), and then eliminate 
events in which the reconstructed mass falls in a narrow window centered on the known top 
quark mass. To accomplish this, we must know the value of the longitudinal momentum for 
the neutrino. We compute this momentum via the on-shell mass relations of the PF-boson 
decay. The quadratic nature of the mass relations produces a twofold ambiguity. In our 
analysis, we include both solutions for the neutrino momentum. After reconstruction, events 
that result in a value of M^i, within a window (of 10, 15, or 20 GeV) around the measured 
top quark mass (we assume rrit = 175 GeV) are rejected. For all three values of the window. 
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FIG. 2: The event rate as a function of ^r for DPS and SPS. On the left, all backgrounds are 
included while the plot on the right compares the DPS events to those from ti alone. While the 
DPS signal is concentrated in the py < 45 GeV range, the majority of the ti background lies above 
this range. Therefore, imposing an maximum JEj^ cut of 45 GeV can greatly reduce the background 
coming from ti production. 

we find a slight improveinent of signal-to-background ratio but accompanied by a decrease 
in the significance {S/yB) of the signal associated with the overall decrease in event rate. 

The poor discriminating power of this mass reconstruction method results from two is- 
sues. First, the twofold ambiguity for the longitudinal component of the neutrino momentum 
provides a combinatorial background. In addition, while the ti events which pass the accep- 
tance cuts are predominantly from the t(th decay mode, about 30% are from the titi mode. 
With two neutrinos in the final state, the on-shell mass relations are not applicable and do 
not provide a unique set of neutrino momenta. For these reasons, the mass reconstruction 
observable is not considered a good discriminator. 

Jets from final states that contain top quarks tend to have a hard spectrum, associated 
with the large top quark mass. A possible observable for tt background rejection is therefore 
the transverse momentum of the leading object (either a jet or a charged lepton). The pt 
spectrum of the leading object tends to be soft in DPS events [5|. However, SPS production 
of Whh yields a rather hard pt spectrum since the bottom quarks recoil against the W^ 
boson. In Fig. [31 we show the pt distributions for the leading object in Whb DPS production 
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and for the remainder (which includes SPS Wbb production). We see that the DPS events 
indeed populate the lower bins of the allowed pr spectrum, while the SPS and background 
events result in a harder spectrum. When we compare the usefulness of placing a cut on the 
upper value of p^ with the improvement we find with the cut on the maximum value of ^t, 
we conclude that the maximum ^t cut provides better significance. If we use both cuts, we 
find that the DPS signal significance is degraded. 



(Teff =12 mb 
L dt = 10 fb"' 




FIG. 3: The event rate as a function of the transverse momentum of the leading object (either a 
jet or lepton). The SPS contribution has a harder px spectrum. 

We conclude this section with the statement that of the three possibilities to reduce the 
ti background we considered, a cut to restrict ^t from above appears to offer the best 
advantage, and it is the only cut we impose in addition to the acceptance cuts specified 
above. 

IV. DISCRIMINATION BETWEEN DPS AND SPS CONTRIBUTIONS TO Wbb 



To separate the DPS events from those of SPS origin, we find it convenient to employ 
quantities which take into account information from the entire final state. One useful ob- 



servable is S'p^, defined as 



i: 




\pT{bi,b2) 



\pTibi)\ + \pTib2)\ 



+ 



\PT{i,^T)\ 



(17) 
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In our case Pt(&1j&2) is the vector sum of the transverse momenta of the two b jets, and 
Pt{^, ^t) is the vector sum of ^t and the transverse momentum of the charged lepton in 
the final state. 
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FIG. 4: The event rate for Wbb production from DPS (left) and SPS (right) as a function of S' . In 
the DPS case, the distribution is peaked toward Sp^ ~ 0; SPS production of Wbb produces bottom 
quarks that are not back-to-back, resulting in a broad distribution and a peak near S' ~ 1. 

In DPS production, the bottom quarks are produced roughly back-to-back such that the 
vector sum of their transverse momenta tends to vanish. Likewise, the vector sum of the 
lepton and neutrino momenta tends to be small (with corrections from the boosted W^). 
Thus, the S' distribution for the DPS process exhibits an enhancement at low S' , as 
shown in Fig. HJ The peak does not occur at exactly S' = owing to NLO real radiation 
that alters the back-to-back nature of the bb and iu systems. On the other hand, SPS 
production of Wbb final states does not favor back-to-back configurations, and it exhibits a 
peak near S' = 1. This feature is linked to the fact that many bb pairs are produced from 
gluon splitting J5[. 

The clean separation in S' between the DPS and SPS Wbb processes exhibited in Fig. H] 
is obscured once the ti background is included (e.g., see the left side of Fig. E]). Figure \5\ 
illustrates the effectiveness of the maximum ^t cut in reducing the ti background in the 
S' distribution. After the cut, a sharp peak is evident in the region of small S' where 
DPS events are expected to reside. 
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FIG. 5: The S' distribution for DPS and SPS production of Wbb including all relevant back- 
grounds. On the left, only the minimal acceptance cuts are imposed, while, on the right, an 
additional maximum ^t cut is imposed {J^t < 45 GeV). Imposing a maximum J<]t cut greatly 
reduces the background and produces a sharp peak in an S' region where DPS is expected to 
dominate. 

The plot on the right side of Fig. El shows that extraction of a relatively clean DPS sample 
can be accomplished by imposing a maximum S' cut. The last column of Tabled shows that 
a cut S'p^ < 0.2 reduces the SPS Wbb rate while leaving the DPS signal relatively unaffected. 
In the end, the major background arises from DPS Wjj, as is expected since this process 
inhabits the same kinematic regions as the DPS Wbb signal. Despite this background, we 
find a statistical significance for the presence of DPS Wbb of S/\fB = 173/a/I97 = 12.3. 

A. Further discrimination 

Observables which take into account the angular distribution of events are also useful in 
the search for DPS. Figure [6] depicts three such observables. In the top-left plot, we show 
the event rates for DPS Wbb and the backgrounds (SPS Wbb included) as a function of the 
angle between the normals to the two planes defined by the bb and iu systems. These planes 
are defined in the partonic center-of-mass frame and are specified by the three-momenta of 
the outgoing jets or leptons. The angle between the two planes defined by the bb and iu 
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FIG. 6: The event rate as a function of the angle between the normals of the two planes defined 
by the bb and iv systems (top- left), the azimuthal angle between the charged lepton and the total 
momentum vector of the bb system (top-right) and the azimuthal angle between the transverse 
momentum vectors of the bb and i ^t systems (bottom). In SPS events, it is apparent that there 
is a strong correlation in the angles. However, there is no such correlation present in the DPS 
events. 

systems is: 

cos /^Qu>,iv = n^ipi, ^2) ■ ^3(^, z/) (18) 

where ns{i,j) is the unit three- vector normal to the plane defined by the i — j system and 
61(62) is the leading (next-to-leading) b jet. In order to construct the normals ^3(61,62) 
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and n^{C,., z/), we require full event reconstruction using the on-shell W^-boson mass relations. 
We see that the distribution of the DPS events is rather flat, aside from the cut-induced 
suppressions at Qi^i^i^, ~ and ~ vr, whereas the SPS events show a strong correlation, with 
a distribution that peaks near A0f,i;£^ ~ |. 

In the top-right plot of Fig. [6l we show the event rates as a function of the azimuthal 
angle between the charged lepton and the total momentum vector of the hh system. No 
information from the neutrino is used. In the bottom plot, we show the event rates as a 
function of the azimuthal angle between the transverse momentum vectors of the hh and 
I ^T systems. Since this azimuthal angle is defined in the transverse plane, it requires only 
^T- Full event reconstruction to determine the neutrino momentum is not needed. In both 
cases, the shape of the DPS distribution is fiat while the SPS distribution shows a strong 
correlation, with a preference for values toward vr. 

In all three plots of Fig. |6l it is clear that DPS and SPS exhibit different behaviors as 
a function of angular observables. However, the dominance of SPS Whh and backgrounds 
over DPS Whh for the full range of these observables makes it impossible to extract a DPS 
Whh signal from these distributions by themselves. 

B. Two-dimensional distributions 

Despite the dominance of the Whh SPS contribution and the backgrounds over the DPS 
Whh contribution, the angular distributions can still be extremely useful when used in con- 
junction with other observables. Two-dimensional distributions of one variable against an- 
other show distinct regions of DPS dominance (or SPS and background dominance). In 
Fig. U\ we construct two such scatter plots. On the left, we show S' versus the angle 
between the normals of the two planes defined by the hh and iu systems (A6(,5£;^), while, 
on the right, we show Sp^ versus the azimuthal angle between the charged lepton and the 
total momentum vector of the hh system. In both plots, we see that the DPS events reside 
predominantly in the lower half of the plane (small S' ) and are distributed evenly in the 
angular variable. The separation between DPS Whh and the SPS component is not as pro- 
nounced in the S'p^ — AO^^^^^ plane as we saw in our earlier study of hbjj [5|. In the Whh 
case, the background events are more evenly distributed over the full plane, to some extent 
resulting from inclusion of both solutions for the neutrino's longitudinal momentum in the 
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FIG. 7: Two-dimensional distributions of events in the variables 5^^ and AGj,^ ^^^ (left) and S'p^ 
and A(j)f^i ^ (right). In both cases, the Whh DPS events (denoted by red x) lie in the lower half of the 
plane, while the Whb SPS and background events (denoted by blue dots) occupy the upper half. 
The plot on the right, in which reconstruction of only the lepton direction is required, appears to 
achieve a cleaner separation, with SPS and background events concentrated in the upper right-hand 
corner of the plane. 

W^ decay. (The greater density of points in the left plot of Fig. [7] relative to the right plot 
is explained by the fact that both solutions for the neutrino momentum are included in the 
left plot). 

As shown in the plot on the right of Fig. [71 the SPS Whb and background events in the 
S' — A(j)),i£ show a strong preference for upper right-hand corner of the plane. This two- 
dimensional distribution indicates that cuts on the 5*^^ and A0^j;^ variables should permit 
extraction of an enriched sample of DPS Wbb events. Inclusion of ^t associated with 
the missing neutrino allows even better separation. In the left plot of Fig. |H1 we show the 
two-dimensional distribution of S'^^ and Acphb^i-^^. This distribution shows a high degree 
of separation between the DPS Wbb and the SPS plus background samples. To quantify 
the degree of separation, we define a region in this plane that gives the highest statistical 
significance. Its boundary is denoted by the black box in the left panel of Fig. |8l Restricting 
^pt ^ ^■'^^ ^^^ ^4>bb,i^T < 37r/4, we find a a sample of 154 signal and 103 background 
events, corresponding to a statistical significance of S/ \fB = 15.2. 
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FIG. 8: The two-dimensional distribution of events in the variables S' and A^^fc^^j, (left). The 
VF66 DPS events are denoted by red x, while the WIA) SPS and background events are denoted by 
blue dots. The box denotes the boundary which gives the highest statistical significance. On the 
right, we show the dependence of the significance as a function of the corners of the box. 

By employing distributions in both S'^^ and t\(^\,h,iYfT-- ^^ achieve a better significance 
than from S' alone. By utilizing only S' we obtain a lower significance of 12.7. In the 
right plot of Fig. [HI we show the dependence of the significance on the placement of the 
box. As long as the maximum value of l^<^hh,e. y/t is in the 7r/2-37r/4 range, a statistically 
significant extraction of DPS Whh from the other events can be obtained, given our assumed 
effective cross section acfr = 12 mb and luminosity. 

We suggest experimental analyses of Whb at the LHC in terms of the two-dimensional 
distributions presented in this section with the goal to establish whether a discernible DPS 
signal is found. Assuming success, the pr dependence of the leading object and other 
properties of these DPS events can be contrasted with those of the remainder to establish 
whether the expected properties of DPS are seen. The enriched DPS event sample can be 
used for a direct measurement of the effective cross section af,<s_. 
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V. CONCLUSIONS AND FURTHER WORK 

In this paper, we investigate the possibihty to observe double parton scattering at the 
early LHC in the pp — ?■ WbbX — )■ iubbX process. Our analysis begins with the basic 
assumption that Wbb production consists of two components: the traditional single par- 
ton scattering process and the double parton scattering process where two individual hard 
scatterings produce the Wbb final state, as depicted in the right panel of Fig. [1] 

After identifying the most relevant background processes, we pinpoint a set of observables 
and cuts which would allow for the best separation between the DPS Wbb signal and the 
backgrounds (including the SPS Wbb process). To provide the most precise predictions 
possible, we generate the DPS Wbb signal event sample, the SPS Wbb sample, and the 
dominant background event samples at next-to-leading order in QCD. The main obstacles 
in the extraction of the DPS signal are the backgrounds from ti production and the SPS 
Wbb component. The most efficient way to suppress the tt background is with an upper cut 
on the missing transverse energy of the event, since top quark decays result in larger values 
of ^T- 

To separate the DPS component of Wbb from the SPS component, we find it useful to 
employ observables which take into account information on the full final state rather than 
observables which involve one or two particles. Examples are the Sp^ variable (defined in 
Eq. ( !T7|) ) and the angle (AG^j^^^) between the two planes defined by the bb and £z/ systems, 
respectively. By displaying the information from these two observables in two-dimensional 
distributions, we show in Sec JIVBl that it is possible to identify distinct regions in phase 
space where the DPS events reside. Utilizing cuts on these observables that enhance the 
DPS Wbb sample, we find that the DPS signal can be observed with a statistical significance 
in the range S/yB ~ 12 — 15. 

A similar study of the DPS and SPS and background contributions to the Zhb final state 
would be a valuable contribution. We remark, however, that the NLO calculation of the 
SPS component of this final state has not yet been implemented in a numerical code such 
as POWHEG needed for a differential analysis like ours for Wbb. 

The focus in the present paper is on establishing double parton scattering as a discernible 
physics process at LHC energies and measuring the size of its contribution. Once DPS 
production of Wbb is observed, it will be interesting to assess its potential significance as 
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a background in searches for other physics, such as Higgs boson production in association 
with a W boson (where the Higgs boson decays as H ^ bb), and precise studies of single 
top quark production where new physics could contribute to the Wtb vertex. A detailed 
analysis of either of these channels would require a different set of optimized physics cuts and 
is beyond the scope of this paper. We limit ourselves here to showing the bb invariant mass 
distribution for the iubb final state in Fig. [91 These results are for illustration only since they 
are based on the cuts outlined in this study. We see that the DPS Wbb component alters 
the overall shape of the bb mass spectrum, enhancing the small mass region. This feature is 
consistent with our earlier observation that the px spectrum of leading jets is softer in the 
DPS component. In Fig. [HI we see that the DPS component contributes primarily in the 
region below 120 GeV or so. At face value, it does not seem to pose a hindrance for searches 
for Higgs bosons in the HW channel. However, Wbb DPS could be a significant background 
in the search for new particles, with masses in the 50 - 100 GeV range and appearing as 
resonances in M^b, and it should be accounted for in any analysis. 

120 r „ 

Wbb DPS 
SPS + bkgs 
Total 
(Teff =12 mb 
Ldt= 10 fb"' 




100 

Mbb(GeV) 



200 



FIG. 9: The event rate as a function of the invariant mass of the bb system using the cuts outlined 
in the text. 
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